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1.  INTRODUCTION 


The  interactions  between  the  exhaust  from  the  shuttle  Primary  Reaction  Control 
System  (PRCS)  engine  firings  and  the  ambient  atmosphere  were  investigated  using 
the  SOCRATES  code.  The  calculations  are  compared  to  ground-based  data 
observed  from  the  Air  Force  Maui  Optical  Station  (AMOS)  [Broadfoot,  et  al. , 
1992].  Data  were  obtained  for  the  shuttle  at  an  altitude  of  320  km,  for  three 
angles  of  attack  (180,  90,  and  0  degrees)  corresponding  to  firings  into  the  ram, 
perpendicular,  and  into  the  wake  of  the  shuttle  direction  of  motion,  respectively. 
The  data  analyses  of  this  experiment  and  some  modeling  efforts  have  been 
investigated  by  Broadfoot,  et  al.  [1992]  and  Elgin,  et  al.  [1991],  The  modeling 
efforts  here  are  focused  mainly  on  simulating  the  general  character  of  the  plumes 
in  the  three  orientations,  and  trying  to  understand  the  characteristics  of  the 
emission  at  6300  A.  This  feature  is  attributed  to  the  atomic  oxygen  red  line 
emission  from  transition  between  the  O^D)  excited  state  and  the  0(3P)  ground 
state.  Two  processes  were  considered  for  formation  of  0(lD):  excitation  of 
atmospheric  0(3P)  by  collision  with  the  exhaust  of  the  shuttle  engines,  and  charge 
exchange  between  O*  and  exhaust  H2O. 

The  SOCRATES  calculations  have  been  done  for  each  of  the  three  angles  of 
attack.  The  following  sections  summarize  this  simulation  work.  Section  2  outlines 
the  SOCRATES  code  algorithms.  In  Section  3,  we  describe  the  AMOS 
measurements  and  the  SOCRATES  model  representation  of  the  experiment. 
Section  4  presents  the  simulation  results  in  the  form  of  temporal  and  spatial 
distributions  of  the  emissions.  Section  5  gives  a  summary  and  conclusion  of  the 
investigation. 


2.  SOCRATES  METHOD  OF  SIMULATION 


The  SOCRATES  code  is  based  on  the  direct  simulation  Monte  Carlo  technique  that 
provides  a  powerful  tool  for  the  simulation  of  real  gas  flows  [Bird,  1981].  The 
SOCRATES  contamination-interaction  program  has  been  developed  to  account  for 
contamination  on  spacecraft  and  to  study  the  neighboring  flow  field  around  the 
shuttle  [Elgin  and  Sundberg,  1988].  The  direct  simulation  Monte  Carlo  method 
has  been  revised  and  extended  significantly  to  account  for  energy  dependent 
collision  cross  sections  and  a  statistical  model  for  internal  energy  effects 
[Borgnakke  and  Larsen ,  1975]  and  [Elgin  and  Sundberg,  1988].  The  method  is 
briefly  described  here.  The  solution  region  is  a  network  of  cells  of  different  sizes. 
Initially,  a  discrete  number  of  molecules  are  stored  in  the  computer,  based  on  their 
velocity  components,  position  coordinates,  and  other  necessary  information.  The 
molecular  motion  and  intermolecular  collisions  in  the  simulated  region  advance  and 
are  modified  with  time,  in  a  two-stage  process.  In  the  first  stage,  the  molecules 
are  propagated  along  their  trajectories  according  to  their  velocity  components  and 
the  time  increment.  At  this  stage,  some  of  the  molecules  may  leave  the  solution 
domain  and  some  will  be  introduced  according  to  the  boundary  conditions  of  the 
problem.  In  the  second  stage,  a  typical  set  of  collisions  will  be  simulated  among 
the  molecules  in  each  cell  with  appropriate  time  increment.  The  repetition  of  these 
two  stages  over  the  small  time  interval  will  uncouple  the  molecular  motion  and 
intermolecular  collisions.  In  a  SOCRATES  simulation,  the  solution  is  periodically 
sampled  by  accumulating  statistical  sums  of  number  densities,  velocities  and  other 
basic  properties.  The  solution  is  run  repeatedly  until  statistical  deviations  are 
reduced  to  a  desired  limit,  and  the  physically  meaningful  output  quantities  are 
computed  from  the  statistical  sums.  The  number  of  molecules  represented  are 
typically  many  thousands  at  a  time,  which  is  much  smaller  than  the  number  of 
molecules  occurring  in  real  flows.  Therefore,  the  construction  of  a  dynamically 
similar  flow  to  be  simulated  in  the  computer  is  an  essential  feature  of  the  method 
[Elgin  and  Bernstein  ,1992]. 

In  the  code,  the  collision  cross  section  is  defined  by  the  variable  hard  sphere  model 
(VHS)  which  depends  on  the  relative  velocity  vr  between  two  molecules.  The 
collision  cross  section  can  be  stated  as 
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where  aref  and  vref  are  the  reference  collision  cross  section  and  velocity, 
respectively  [Elgin,  et  al,  1990],  oj  is  a  constant  parameter  that  has  a  value  of 
0.25  for  this  investigation. 

Reactive  collisions  between  molecules  and/or  atoms  can  be  simulated  directly  by 
SOCRATES.  The  Monte  Carlo  program  simulates  the  reaction  with  a  probability 
that  is  related  to  the  ratio  of  the  reactive  to  collision  cross  section  at  the  relative 
velocity  for  the  collision  [Elgin  and  Sundberg ,  1992]. 

The  present  work  employs  the  option  of  using  the  Arrhenius  rate  constant  in  the 
code  to  calculate  the  reactive  collision  cross  section.  The  rate  constant  has  the 
form  of 


kr  =  a  r  *  R°T 


(2) 


where  A  and  n  are  constant  parameters,  R0  and  T  are  gas  constant  and 
temperature,  respectively,  and  Ea  is  the  activation  energy.  The  reactive  cross 
section,  a *,  corresponding  to  Equation  (2),  depends  on  the  relative  collision  energy 
and  is  given  by  [Bird,  1981] 
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Here,  Ec  is  the  collision  energy,  6y  is  unity  for  like  reactants  and  zero  for  unlike 
reactants,  and  F  is  the  gamma  function.  In  the  case  of  collision  between  two 
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molecules,  the  reactive  cross  section  is  calculated,  and  the  reaction  is  counted  with 
a  weighting  factor  Wr.  The  weighting  factor  is  given  by 

W  =  W— -  (4) 

r  c  a 


where  a *,  a,  and  Wc  are  the  reactive  cross  section,  collision  cross  section,  and 
collision  weighting  factor,  respectively.  Further  discussion  of  the  reactive 
collisions  can  be  found  in  the  report  by  Elgin  and  Sundberg  [1992]. 


3.  DESCRIPTION  OF  THE  EXPERIMENT  AND  CONFIGURATION 
MODEL 


The  experiment  was  based  on  firing  the  Primary  Reaction  Control  System  engines 
of  the  space  shuttle.  The  PRCS  liquid-fuel  engines  were  fired  in  the  ram,  wake, 
and  perpendicular  directions  (angles  of  attack  of  180°,  0°,  and  90°  respectively). 
The  experiment  was  performed  on  October  7,  1990  at  1326  UT,  while  the  Shuttle 
overflew  the  Air  Force  Maui  Optical  Station.  The  shuttle  was  in  the  300-km 
altitude  range.  The  imagery  data  from  PRCS  engine  firings  were  obtained  using 
an  Intensified  Charged-Coupled  Device  (ICCD)  spectrograph  at  the  focus  of  a  200- 
mm  aperture,  f/5  telescope  mounted  on  and  coaligned  with  a  1.6  m  telescope.  The 
details  and  discussion  of  data  acquisition  and  analysis  of  the  experiment  can  be 
fouod  in  reference  [Brocujfoat,  et  al.,  1992]. 

The  main  purpose  of  our  investigation  is  to  simulate  this  experiment  using  the 
SOCRATES  code.  Therefore,  considerable  care  has  been  given  to  use  a  realistic 
set  of  species,  chemical  and  photochemical  reactions,  and  parameters  for  the  PRCS 
<*"gine  and  ambient  atmosphere.  The  solution  domain  for  the  simulations  is  taken 
to  be  a  space  of  dimensions  21  by  20  by  20  km  consisting  of  3332  unevenly 
spaced  cells.  In  each  cell  there  are  3  to  20  molecules  (depending  on  density  of 
species  and  ceil  location)  for  each  species.  This  setup  has  generally  been  retained 
for  all  three  firings  in  the  ram,  perpendicular,  and  wake  directions.  All  of  the 
species  and  their  reactions  and  most  of  the  important  parameters  which  have  been 
used  for  the  simulation  are  given  in  Tables  1-4. 
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The  exhaust  and  atmospheric  mole  fractions  of  the  species  used  in  the  calculations 
are  given  in  the  Table  1.  Since  the  other  species  present  in  the  atmospheric 
composition  have  no  significant  effect  on  the  problem,  they  have  not  been 
considered  in  the  calculations. 


TABLE  1.  Concentrations  of  Species  in  the  SOCRATES  Simulation 


Species 

Atmospheric  Mole 
Fraction 

Exhaust  Mole 
Fraction 

0 

0.9076 

0.0 

N2 

0.0919 

0.309 

0* 

0.0005 

0.0 

h20 

0.0 

0.332 

h2 

0.0 

0.187 

CO 

0.0 

0.136 

co2 

0.0 

0.036 

Table  2  shows  a  list  of  the  species  used  in  the  simulations.  These  species 
represent  engine  exhaust,  as  well  as  atmospheric  ambient  compositions  and 
reactions  products.  The  table  also  lists  the  reference  collision  cross-section  aref, 
between  two  members  of  this  species,  corresponding  reference  relative  collision 
velocity  vref,  number  of  internal  degrees  of  freedom  vif  and  heat  of  formation  for 
the  species. 
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TABLE  2.  Molecular  Parameters  Used  in  the  Calculations 


Species 

aref 

(cm5) 

vref 

(cm/s) 

Heat  of 
Formation 
(kcal/mole) 

vi 

0 

1.750  x  10-13 

2.490  x  105 

59.60 

0.00 

3.540  X  10"15 

2.130  X  105 

0.00 

2.00 

o* 

1.395  x  10-14 

2.490  X  103 

372.40 

0.00 

h20 

3.270  x  10-'5 

2.290  x  103 

-57.80 

3.76 

h2 

1.860  x  10'15 

7.740  x  103 

0.00 

3.29 

CO 

3.460  x  lO'13 

2.090  x  103 

-26.42 

2.00 

co2 

4.330  x  10‘13 

1.710  x  105 

-94.10 

3.58 

0(lD) 

1.750  x  10-'5 

2.490  x  103 

104.90 

0.00 

OH 

1.770  x  10'15 

8.000  >  105 

9.40 

2.00 

h2o+ 

3.270  X  10'3 

2.290  X  103 

232.70 

3.76 

Table  3  lists  the  chemical  reactions  used  in  the  simulations.  There  are  generally 
two  types  of  mechanifms  for  generating  O^D)  given  in  the  table: 


the  charge  exchange  reaction 


h2o  +cr  -  h2o+  +  cx}d) 


and 

0?P)  +  X  0(lD)  +  X 


where  X  is  any  exhaust  species  with  sufficient  kinetic  energy  to  excite  the  0(1D). 
For  example,  in  Table  3,  X  represents  N2,  0(3P),  and  H20.  The  O^D)  can 
undergo  four  quenching  reactions  as  shown  in  Table  3,  which  do  not  produce 
photons.  The  photochemical  decay, 
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0('D)  -  0(3P)  +  hv  (X  =  6300  A) 


is  proposed  to  generate  the  photon  emission  at  6300  A.  The  table  also  shows  the 
values  used  for  the  Arrhenius  rate  coefficients  for  the  different  reactions.  The 
discussion  of  the  constants  in  Table  3  can  be  found  in  the  reference  [Broadfoot, 
etal .,  1992]. 


TABLE  3.  Rate  Coefficient  Used  in  SOCRATES,  k=  A  Tn  exp(-Ea/kT) 


Ea 

Reaction 

(kcal/ 

AH 

Reaction 

A 

n 

mole) 

(kcal/mole) 

0(3P)  +  N2  -»  0(‘D)  +  N, 

or3P)+o(3P)  -»  o(‘d)+o(3p) 

1.28  x  10-“ 

0.0 

45.3 

45.3 

4.9  x  10-‘2 

0.0 

45.3 

45.3 

0(3P)  +  H20  -  O(’D)  +  HjO 

1.28  x  lO'12 

0.0 

45.3 

45.3 

0(‘D)  +  Nj  -  0(3P)  +  N, 
0('D)+0(3P)  -»  0(3P)+0(3P) 

2.3  x  10-“ 

0.0 

0.0 

-45.3 

8.0  x  10’12 

0.0 

0.0 

-45.3 

0(‘D)  +  HjO  -*  0(3P)  +  HjO 

2.3  x  10'“ 

0.0 

0.0 

-45.3 

O('D)  +  HjO  -»  OH  +  OH 

2.2  x  lO'10 

0.0 

0.0 

-28.3 

HjO  +  0+  -»  H20+  +  O('D) 

3.23  x  lO’9 

-0.03 

23.0 

23.0 

0(!D)  •+  0(3P)  (photochemical) 

1.5  X  102  (lifetime, 

seconds) 

Table  4  represents  the  summary  of  information  for  the  ambient  environment  and 
PRCS  engine  exhaust  bums  released  into  space.  The  first  four  parameters  given 
in  the  table  show  the  ambient  conditions,  and  the  rest  are  exhaust  system 
parameters. 
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TABLE  4.  Parameters  Used  in  the  Simulations 


Altitude 

Velocity  of  Spacecraft 

Atmospheric  Temperature 

Atmospheric  Number  density 

Mass  Flow 

Exit  Plane  Area 

Ratio  of  Specific  Heat 

Exit  Mach  Number 

Exit  Nozzle  Half  Angle 

Exit  Plane  Density 

Exit  Plane  Number  Density 

Exit  Plane  Velocity 

Exit  Plane  Speed  of  Sound 

Exit  Plane  Temperature 

Exit  Plane  Pressure 

Thrust 


320  km 
7.4  km/s 
925°  K 

5.74  x  108  cm'3 
1.227  X  103  g/s 
4.66  x  102  cm2 
1.297 
4.02 
8.0° 

8.872  x  10'6  g/cm3 
2.618  x  1017  molecules/cm3 
2.969  X  105  cm/s 
7.386  x  104  cm/s 
1.032  x  1030  K 
3.731  x  104  dynes/cm2 
8.584  x  102  lb 


4.  DISCUSSION  OF  RESULTS 


The  results  discussed  here  are  based  on  a  time  dependent  simulation  using  the 
SOCRATES  code  for  calculating  the  distribution  of  0(1D)  -»  0(3P)  photons.  The 
problem  was  simulated  in  20  time  steps,  spanning  8.0  seconds  beginning  with  the 
start  of  the  engine.  It  was  assumed  that  the  PROS  engine  was  shut  down  at  3.0 
seconds  after  it  started.  The  plots  presented  here  are  a  representative  selection  of 
the  20  step  times.  The  simulations  were  performed  in  three  directions  of  ram, 
perpendicular,  and  wake. 

Figures  1-3  show  the  contour  plots  for  the  evolution  of  the  6300  A  plume  in 
radiance  (Watts/steradian/cm2)  for  0(lD)  0(3P)  transition  as  a  function  of  time 

for  ram,  perpendicular,  and  wake  burns,  respectively.  In  these  figures,  the 
calculated  results  are  plotted  for  six  step  times:  two  contour  plots  before  shutdown 
of  the  engine,  one  at  the  time  of  shutdown  (3.0  see),  and  three  contour  plots  after 
the  shutdown  of  the  engines.  As  can  be  seen  from  the  figures  (especially  the  ram 
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Figure  1.  Contour  plots  showing  the  calculated  6300  A  plume  radiance  for  a  ram  burn  at  several  time  steps. 
The  shuttle  is  located  at  the  origin  of  coordinates,  and  the  wind  is  flowing  to  the  right. 
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Figure  2.  Contour  plots  showing  the  calculated  6300  A  plume  radiance  for  a  perpendicular  burn  at  several  time 
steps.  The  shuttle  is  located  at  the  origin  of  coordinates,  and  the  wind  is  flowing  to  the  right. 
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Figure  3.  Contour  plots  showing  the  calculated  6300  A  plume  radiance  for  a  wake  bum  at  several  time  steps 
The  shuttle  is  located  at  the  origin  of  coordinates,  and  the  wind  is  flowing  to  the  right. 


case,  Figure  1),  the  relatively  slow  radiance  decay  after  the  engine  shutdown  is 
due  to  the  long  radiative  lifetime  of  O^D).  The  emission  disappears  mostly  due 
to  diffusion  of  O^D). 

Figures  4-5  show  the  radiant  intensity  (Watts/steradian)  for  6300  A  emission  as  a 
function  of  time  in  the  three  directions  of  ram,  perpendicular,  and  wake.  Figure 
4  is  plotted  based  on  the  SOCRATES  calculation  of  radiant  intensity  in  Watts/(sr 
cm2).  Using  an  auxiliary  program  called  MAP,  this  data  was  integrated  through 
the  line  of  sight  to  give  radiance  in  Watts/sr.  The  radiant  intensity  in  Watts/sr  can 
also  be  calculated  from  the  SOCRATES  output  of  the  O^D)  number  density. 
Using  the  MAP  program,  the  total  number  of  excited  molecules,  N,  in  the  field 
of  view  is  obtained.  From  this,  we  can  calculate  the  average  decay  rate  as  N/r 
where  t  is  the  lifetime  for  the  excited  state.  The  energy  emitted  by  each  photon 
is  Ae  =  hc/X  where  h  =  Plank’s  constant,  c  is  the  speed  of  light,  and  X  is  the 
wavelength.  Hence,  the  total  energy  emitted  per  second  is  e  =  N/r  Ae  and  the 
energy  per  steradian  is  E=  e  /  (4  it).  In  Watts/sr,  E  -  (  N/r  1/X  1.99  x  10'19), 
where  X  is  in  micron.  Figure  5  is  based  on  the  density  variations  of  0(*D)  vs. 
time  in  the  directions  of  ram,  perpendicular,  and  wake.  The  calculation  proceeded 
as  described  above  using  the  density  and  life  time  of  0(lD),  and  a  converting 
factor  to  obtain  the  radiant  intensity  (Watts/steradian).  Figures  4  and  5  are  almost 
identical,  as  expected.  Figure  6  shows  the  decay  rate  (reactions/sec)  obtained  from 
the  MAP  program  for  the  reaction  0(lD)  -*  0(3P)  photons.  This  plot  is 
proportional  to  Figure  4,  but  given  in  different  units.  It  should  also  be  noted  that 
Figures  4-6  were  calculated  using  the  field  of  view  of  20  by  20  km. 

Figures  7-12  show  the  same  cases  as  for  Figures  1-3,  but  depicted  in  color  to 
better  represent  the  morphology  of  the  plumes.  The  color  plots  were  created  using 
the  IRMA  program  [Tautz,  1993].  This  code  can  read  the  SOCRATES  output 
files  to  extract  the  reaction  data,  and  then  integrate  along  an  arbitrary  line  of  sight 
to  compute  the  plume  radiance.  There  are  two  sets  of  plots  for  each  case  (ram, 
perpendicular,  and  wake  firing).  The  first  set  of  four  frames  shows  the 
development  of  the  plume  up  to  the  engine  shut  off  time  (3  sec.)  and  the  second 
set  shows  the  dispersion  of  the  plume  after  shut  off.  In  each  plot,  the  shuttle  is 
located  at  the  origin  of  the  (x,y)  coordinates  indicated  in  the  frames.  The  viewing 
window  is  16  by  16  km,  and  the  length  of  the  x-y  axis  is  1  km.  In  the  simulation, 
the  shuttle  velocity  is  in  the  negative  x  direction,  so  that  in  the  shuttle  rest  frame 
the  ambient  environment  flows  from  left  to  right.  A  common  color  legend, 
representing  radiance  in  Watts/(sr  cm2),  has  been  used  for  all  the  frames.  It  is 
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Figure  5.  This  figure  is  seme  as  Figure  4  except  that,  it  is  calculated  from  the  number  density  and 
lifetime  of  0(*D). 
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Figure  7.  Photon  emission  at  6300  A  for  ram  burn  before  shutdown- 


Photon  emission  at  6300  A  for  rant  burn  after  shutdown. 
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Figure  9.  Photon  emission  at  6300  A  for  perpendicular  burn  before  shutdown. 
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Figure  10.  Photon  emission  at  6300  A  for  perpendicular  burn  after  shutdown. 
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Figure  11.  Photon  emission  at  630*3  A  for  wake  bum  before  shutdown. 


Figure  12.  Photon  emission  at  6300  A  for  wake  bum  after  shutdown. 


evident  that  the  ram  burn  produces  the  most  intense  plume,  followed  by  the 
perpendicular  burn,  and  that  the  wake  burn  is  barely  visible  on  this  scale. 

The  results  that  have  been  discussed  thus  far  are  concerned  with  the  O^D)  photon 
emission.  However,  for  better  presentation  of  the  problem,  it  is  important  to  study 
the  behavior  of  the  other  species  and  their  reactions  in  the  simulation.  Therefore, 
a  brief  discussion  of  the  simulation  for  the  rest  of  the  species  which  are  illustrated 
in  Tables  2  and  3  will  be  given  here.  Figures  13-22  show  an  arbitrary  z-plane 
slice,  near  the  center  of  the  interaction  region  (z  =  0.259  km)  of  species  densities 
and  reactions  rates  from  the  SOCRATES  output  for  the  ram  case  at  2,8  sec  after 
the  engine  has  been  fired.  The  figures  are  plotted  in  the  field  of  view  of  20  by  20 
km. 

Figures  13-18  show  the  number  densities  of  the  atmospheric  and  engine  species, 
as  well  as  the  species  which  are  produced  by  the  reactions  shown  in  Table  3.  In 
the  figures,  four  features  can  be  seen.  First,  the  atmospheric  species  (0,  0+,  and 
N2)  are  depleted  starting  from  the  source  (the  shuttle’s  engine)  and  extending 
toward  the  wake  of  the  shuttle  for  more  than  5  km.  For  example,  Figure  13  (a) 
shows  the  density  distributions  of  atmospheric  O  atoms.  It  is  evident  that  the  O 
density  has  been  depleted  at  the  location  of  the  shuttle  engine  and  extends  toward 
the  wake  of  the  shuttle.  This  depletion  is  produced  by  a  combination  of  the 
chemical  reactions  in  the  ram  and  deflection  of  the  O  atoms.  Second,  the  density 
of  0  atoms  is  increased  in  the  ram  direction  due  to  the  ‘snow  plow’  effect.  From 
Figure  13  (a),  it  can  be  seen  that  the  density  is  enhanced  about  15%  above  ambient 
values.  Third,  the  exhaust  species  are  concentrated  mainly  on  the  ram  side  of  the 
shuttle.  Fourth,  N2  is  a  common  species  between  the  atmosphere  (9.2%)  and 
exhaust  (30.9%)  compositions,  and  Figures  15  and  16  show  that  the  N2  is  acting 
both  like  exhaust  engine  species,  as  well  as  an  atmospheric  species,  as  it  shows 
some  downstream  depletion.  In  the  same  figures,  the  H20,  which  only  exists  in 
the  exhaust  composition,  does  not  show  any  depletion.  It  should  be  noted  that 
Figures  15  (a)  and  (b)  are  plotted  based  on  the  maximum  and  upper  values  of  the 
densities  and  Figures  16  (a)  and  (b )  are  plotted  based  on  a  lower  density  range. 
As  discussed  earlier,  the  spatial  distribution  of  the  reaction  rate  for  O^D)  -*  0(*P) 
emissions  can  be  estimated  by  dividing  the  O^D)  density,  given  in  Figure  13  (b), 
by  the  lifetime  of  150  seconds. 
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0+  density  (atoms/cm*3)  at  2.8  seconds  H20+  density  (molecules/cm*3)  at  2.8  seconds 
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Figure  14(a).  Charged  O  densities  on  an  arbitrai-y  Figure  14(b).  Charged  H20  densities  on  an 

central  z  plane  slice.  arbitrary  central  z  plane  slice. 
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Figure  15(a)*  N2  and  densities  on  an  arbitrary  Figure  15(b).  H20  densities  on  an  arbitrary  central 

central  z  plane  slice,  plotted  based  on  the  z  plane  slice,  plotted  based  on  the  maximum  and 

maximum  and  upper  values  of  the  densities.  upper  values  of  the  densities. 
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C02  density  (moleaj!es/cm*3)  at  2.8  seconds  OH  density  (mo!ecuIes/cm*3)  at  2.8  seconds 
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Figure  18(a).  C02  and  densities  on  an  arbitraiy  Figure  18(b).  OH  densities  on  an  arbitrary  central 

central  z  plane  slice.  z  plane  slice. 
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Figure  19.  Reaction  rates  on  an  arbitrary  central  z  plane  slice,  for  reaction  given  above. 


Figure  22.  Reactions  rates  on  an  arbitrary  central  z  plane  slice,  for  reaction  given  above. 


Figures  19-22  show  the  reaction  rates  (#/sec  cm3)  for  the  reactions  in  Table  3. 
Figures  19  and  20  give  the  spatial  distribution  for  the  reactions  that  generate 
excited  O^D)  atoms.  The  figures  show  that  the  dominant  reaction  rates  occur  in 
Figures  19  (a)  and  20  (a),  where  an  0  atom  collides  with  an  N2  molecule  or  a 
water  molecule.  Figures  21  and  22  represent  the  spatial  distribution  for  the  four 
quenching  reactions  given  in  Table  3.  It  can  be  seen  that  Figure  22  (b),  where  an 
0  atom  collides  with  a  water  molecule,  shows  the  most  intense  quenching  rate. 
The  reaction  rates  generating  O^D)  atoms  are  stronger  than  the  quenching  rates, 
so  that  the  O^D)  density  is  increasing  at  time  2.8  seconds.  The  most  intense 
generating  reactions  are  concentrated  near  the  shuttle  and  the  excited  0(JD)  atoms 
diffuse  out  to  form  the  plume. 


5.  SUMMARY  AND  CONCLUSION 


The  SOCRATES  code  was  used  to  simulate  the  interactions  between  the  PRCS 
engine  firings  and  the  ambient  atmosphere  in  the  three  orientations  of  ram, 
perpendicular,  and  wake  of  the  shuttle  motion.  The  calculations  were  made  at  an 
altitude  of  320  km.  The  simulation  results  indicated  that  when  the  exhaust  of  the 
space  shuttle  interacts  with  the  ambient  atmosphere,  it  generates  long-lasting 
emissions  at  6300  A.  The  intensity  of  these  emissions  varied  from  strong  to  weak 
for  ram  to  wake  firings,  respectively.  The  discussions  of  the  intensity  of  0(*D  -* 
3P)  transition  photons  as  a  function  of  firing  directions  (ram,  perpendicular,  and 
wake)  and  time  were  illustrated  in  Figures  1  through  6.  A  series  of  color  plots 
was  presented  to  show  the  plume  developments  as  a  function  of  time  for  the  three 
orientations  of  engine  firing  (Figures  7  through  12).  A  series  of  figures  was  also 
shown  to  illustrate  the  involvement  of  other  species  and  their  reactions  in  the 
simulation  (Figures  13-22). 
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